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The problem of wave propagation in stratified media is reexamined from the stochastic 
point of view by introducing the concept of random media. For a small inhomogeneity, 
this problem is investigated by utilizing the effective refractive coefficient of Chen [1964]. 
For a large inhomogeneity, this problem is treated by subdividing the medium into many 
parallel homogeneous layers of random media and utilizing the derived reflection and trans- 
mission coefficients of Chen [1964]. In each case, it is found that the small randomness 
has no drastic and unexpected effects on the behavior of the wave motion in the media 
considered. Hence, our reliance on the results obtained by a nonrandom approach is 
increased. 

1. Introduction 

The subject of wave propagation through a nonrandom medium has been investigated 
thoroughly by many authors [Brekhovskikh, 1900; Wait, 1902]. Since, in practice, the 
properties of the medium are either impossible to be measured accurately or subjected to 
random variations due to noise, humidity, wind velocity, thermal fluctuations, etc., this 
problem should be reexamined from the stochastic point of view. Hence, the concept of 
random medium is introduced. 

The wave motion in a transmission medium is described by a vector-valued function 
u(r, t) of the position vector r and the time t. As a consequence of the physical laws govern- 
ing the wave motion, u(r, t) satisfies certain partial differential equations of symmetric hy- 
perbolic type. The effect of the transmission medium on the wave motion is characterized 
by a vector-valued function N(r, t) the "propagation coefficient. " For a random medium 
the coefficient N(r, t, q) depends also upon a parameter g, ranging over a space 12 in which 
a probability density P(q) is defined. If u(f, t, q) describes the wave motion in a random 

medium, then the mean value of u(r, t, q) is defined by (u)= P(q)u(7, t, q)dq. 

In order to distinguish the many phases of random media, we write N(r, t, q)=k(r, t) 
n (r, t, q). Then "homogeneous continuous random medium" is defined as k(r, t) being a 
constant. "Inhomogeneous continuous random medium" is defined as &(r, t) being a con- 
tinuous function of r and t. Finally, "discontinuous random medium" is defined as k(r, t) 
being a discontinuous function of r and t, and n(r, t, q) may have different P(q)'s on the two 
sides of the discontinuity. By the above definitions, we observe that the random stratified 
media can be either inhomogeneous continuous random media or discontinuous random media. 

Many works have been done on the wave propagation in homogeneous continuous random 
media [Chernov, 1960; Keller, 1960 and 1963; Furutsu, 1963]. However, their results fail to 
apply to the problem of wave propagation in random stratified media. It has not been until 
recently in a paper by Y. M. Chen [1964], that the effective refractive coefficient of a slightly 
inhomogeneous continuous random medium and the reflection and transmission coefficients 
for the wave propagation in a discontinuous random medium have been obtained; now one 
can honestly investigate the problem of waA^e propagation in a random stratified medium. 

In this paper we shall study the problem of wave propagation in stratified media of small 
randomness. For a small inhomogeneity, this problem can be treated by utilizing the effective 
refractive coefficient of Chen [1964]. The solutions of this type are obtained for two well- 
known special profiles of the propagation coefficient [Landau and Lifshitz, 1958]. For a large 
inhomogeneity, this problem can be effectively treated by subdividing the medium into a num- 
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ber of parallel homogeneous layers of random media and utilizing the derived reflection and 
transmission coefficients of Chen [1964]. In fact, by taking a sufficiently large number of 
such layers of decreasing thickness, any desired degree of precision may be obtained. In each 
case, the results are compared with those for the nonrandom media and no significant differences 
are found. For simplicity, only scalar waves of harmonic time dependence, e~ iut , are 
considered. 

2. Reflection and Transmission of Plane Waves From Horizontally Stratified 
Random Media of Small Inhomogeneity With Special Profiles 

If r denotes a point in three dimensional space, then u(r, q), characterizing the motion of 
plane in a random stratified medium, satisfies the following reduced wave equation 

vMr,q)+k\(l + ef(T))(l + ew(r,q))] 2 u(r,g)=0 (1) 

where e is a small parameter, /(r) is a continuous function of r and w(r, q) is a continuous function 
of r and q. From Chen [1964], by assuming (w)=0, it is found that up to and including 
terms of order e 2 , the mean value of u(r) satisfies the following differential equation: 

V 2 (u)+k 2 (l + ef(?)y(u)=0, (2) 

k 2 =k 2 (l + e 2 (w 2 )K), (3) 

where 

K=l-2ik r {e mr -l)C{r)dr, (4) 

and the correlation coefficient 

CQr-r'\)=^p (5) 

is assumed to be a function of the distance \r— r'\ only. 

Even for k being real and positive, the imaginary of k can be shown to be positive and its 
real part can be shown to be greater than (l-\-e 2 (w 2 )E?)k [Keller, 1963]. Hence the amplitude 
and phase velocity of the coherent wave, (u), is exponentially attenuated and reduced re- 
spectively by the randomness of the medium. 

Now, f(r) is assumed to be a function of y only, and u(r)=u(y) represents a plane wave, 
coming from — °°, propagating along the ?/-axis. Hence, upon omitting the term e 2 J(y) 2 , (2) 
becomes 

^(u)+Hl+2ef(y))(u)=0. (6) 

If (l+2e/(y)) approaches constants U 2 and V 2 as y-^—co and ?/^+oo respectively (U 2 may or 
may not equal V 2 ), then 

{u)~Ce ilvv asy-xx), (7) 

and 

(u)^Ae ikUv +Be- tkUv as y->— oo . (8) 

In this case, the reflection and transmission coefficients are defined as 

R=A~ 1 B (9) 

and 

T=C (10) 

respectively. 

Case (a): ( 1 + 2 ^(2/))=l- 1 "3^- (11) 
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In order to find the reflection coefficient, one has to find a solution of (6) which has the 
form 

(u) ^constant • e d(l -* )h as y-*x> . (12) 

By introducing a new variable 

p=-e~ aV (13) 

and seeking a solution of the form 

(u)=p- da -* ) *«- 1 *(p) (14) 

(where ^(p) tends to a constant as p—>0, i.e., y-^ 00 ), one finds that *&(p) satisfies the following 
hypergeometric differential equation 

p(l-^)^ ,/ + [l-2^(l-e) 1/2 a- 1 ](l-^)^ / -€Fa- 2 ^ = 0, (15) 

which has as its solution the hypergeometric function 

¥=AfF(*[l — (1— e) va ]&^ (16) 

where M is an arbitrary constant. This function satisfies the imposed condition, i.e., as 
p->0, ^— >1. By using the asymptotic form of ^ as p^> — oo, one obtains the asymptotic 
form of (u) as y— >— <*> 

(u)^(-l)- rkil -^ a ~ l [D 1 e ury +D 2 e- ik \ (17) 

where 

D __ T(-2ika-i)T(-2ika-i(l-e)$+l) 

r[-^(l + (l-e)^- 1 ]r[-^(l + (l-e)^)a- 1 + l] 
and 

D T(2ika- l )T(-2ika- l g-e^+l) 

r[^(i-(i-€)^- 1 ]r[^(i-(i-e)^a- 1 +i] 

One*finally obtains 

^_A^ r(2^a- 1 )r[-^(i+(i-6)^- 1 ]r[-^(i + (i-e)^- 1 +i] 
A r(-2ifo- 1 )r[^(i-(i-e)^- 1 ]r[^(i-(i-e)^)a- 1 +i] 

Since Im (k— fc)>0 and Re &>(l + € 2 (w 2 )iT)Re k [Keller, 1963], if one keeps away from the 
poles and zeros of the gamma functions, the above reflection coefficient R is just an analytic 
continuation of the reflection coefficient for this medium with randomness removed. For all 
the practical cases, the arguments of the above gamma functions do not equal °°, 0, —1, —2, 
— 3, . . . . Thus, the small randomness of the medium only changes and shifts the amplitude 
and phase of the reflection coefficient respectively. 

Case(b): (1+2^))^!- ^/^,^ (21) 

In order to find the transmission coefficient, one has to find the solution of (6). By 
making the substitution 

(tt) = [cosh (2//3- 1 )]" 2X ^(p), (22) 

where 

\=![(l_4/3 2 Fe) 1/2 -l] (23) 

and 

^=-sinh 2 (ypr 1 ), (24) 

one finds that V(p) satisfies the following hypergeometric differential equation 

p(l-p)*' / + [±-(l-2\)p]*'-(\ 2 +l(3 2 k 2 )*=0. (25) 
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Thus 

<«> = Ci[cosh (ytr^-^-l+ihW-l-iSkfrb-smh* (yfi- 1 )] 

+ a[cosh (yp- 1 )]- 2 * sinh (yp-^Fl-^+fktf+l-X-ifkp+h 1|; -sinh 2 (y/T 1 )]. (26) 

The coefficients d and C 2 are determined from the condition that (u)~e m as y— >+co. 
Upon using the proper asymptotic forms of the hypergeometric functions, one obtains 

{u)^(-l)^l(C 1 J 1 -C 2 J 2 )(-h)-^e i ~''+(C l K 1 -C2K 2 )(-l)^e-^] as y-*-», (27) 
and 



where 



(«)~[(0 , 1 «/i+C 2 J 2 )(J)-*e-^+((7 1 Z- 1 +Cyir !! )(|)^ l£ ''] as 2/^+°°, (28) 

J X = =**=*> , (29) 

J 2 = ^=^ =- (30) 

2T(-\+i-ilkp)T(\+l-iik0) 

*« !^* ^, (31) 

r(-x+^^)r(x+i+iifc/5) 

K^ ^* (32) 

2T(-\+i+ikifi)T(\+l+iikp) 

The condition that as y^> + <*< only the transmitted wave is present leads to the following 
relation between (7i and C 2 

Ci«7,+ ftj«=0. (33) 

Finally, 

(C^+CiKi) (-i)- 2 V IV* 3 



and 



T=- 



(4) 



\C\J\ C2J2) 

' KJ 2 -J,K 2 






Zd \d 2 



(34) 



As in case (a), if one keeps away from the poles and zeros of the gamma functions, the above 
transmission coefficient T is just an analytic continuation of the transmission coefficient for 
this medium with the randomness removed. Hence, the small randomness of the medium 
only changes the amplitude and shifts the phase of the transmission coefficient. 

3. Reflection of a Plane Wave From Horizontally Stratified Random Media 

of Large Inhomogeneity 

Any stratified random medium of large inhomogeneity can be approximated by a sufficiently 
large number of parallel homogeneous layers of random media to any degree of accuracy. 
First, the same problem for the nonrandom media is formulated. For the time being, the 
number of layers is taken to be M and they lie horizontally below the x—z plane (fig. 1). A 
plane wave, 

^. =e iJcQX sin 0-f ik y cos C%%\ 

is incident at an angle 6 on the upper surface of the first layer. The wave motion in the mth 
layer below the x—z plane is characterized by the solution of one of the following partial dif- 
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ferential equations, 
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Figuhe 1. 77? e structure of a stratified random 
medium consists of M layers of homogeneous random 
media. 



(V 2 +kl)u n =0, m = l, 2, 3, 4, . . .,M, 



(36) 



where k m is the propagation constant with nonnegative imaginary part, and the boundary 
conditions at the interfaces y=0, yi, y 2 , y?>, • . ., Vm-\ (fig- 1) are expressed by 






dy 



l^my 



&ty=y n -i m=l,2, 3, 



,M, 



(37) 



where f$ m is a constant determined solely by the properties of the rath layer. 

Because of the randomness in these media, it is more convenient to discuss (u m ), the mean 
of u m than u m itself. From Chen [1964], if one assumes (w m )=0 and k m -i^k m , ra=l, 2, 3, . . ., 
M (this implies that the division of medium has to be very fine such that (R 2 ) and (T 2 ) of (85) 
of Chen [1964] are negligible, one finds that up to and including terms of order e 2 (perturbation 
parameter), away from the boundary surfaces (u m ) satisfies (2), (3), (4), and (5) with the sub- 
script m inserted for everything except e, r, and r, and at the boundary surfaces, (u m ) obeys (37) . 

Now we have the form of the general solution as 



(uJ^Ome^+ime-^e**, m=l, 2, 3, . . .,M, 



(38) 



where km=s*-\-^2, and X can take any value. Upon imposing the outgoing wave condition on 
(38) for m=M and inserting (35) into (38) for ra=0, it follows that b M =0 and a =l, s = 
k cos 6, \=k Q sin 6. Upon applying the boundary conditions (37) to (38), the coefficient b is 
obtained as 

, _ — ikpfio cos 6— z\ 

— ifcoftjCOSfl + Z! 
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where 



Z2— Sift tan Si/^i 
Sift +2 2 tan s^i 



Z 2 = S 2 P2 



Z3—S2P2 tan 5 2 ^2 
$2^2+23 tan S 2^2 



^ £?w+i ^mPm tan s m /k m 

o m Pm~r Zm+i ^an S m A6 m 



2ikf- 1 — — ISm- iPm- 1 



SmPmt %-ifcf-i tan s M -\I^m-i 
Sm-iPm-i\SmPm tan SM-\h>M-i 



(39) 



and A w is the thickness of the mth layer. 

Since \s m — s m \ is very small, if one keeps away from the singularities of tan (s m h m ) and of 
the denominator of z m one would expect no drastic change in b . It is also interesting to notice 
that in this case b is not an analytic continuation of the reflection coefficient for this medium 
with randomness removed as in section 2. 
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